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ABSTRACT 

This work analyses the properties of groups of galaxies in the surroundings of clusters. 
On the basis of a very large public Virgo Consortium Simulation, we identified systems 
of galaxies in a wide range of masses. Systems with masses greater than Mcut = 
4 X 10^"* Mq h^^ are considered "host", whereas smaller systems are taken as groups. 
Our results show that groups properties are affected by the proximity of massive hosts. 
Physical properties such as velocity dispersion, internal energy (E) and virial radius, 
show an increment, whereas the mean density decreases as the host-group distance 
is smaller. By analysing groups with different properties, we find that the low mass 
and the weakly bounded {E > 0) subsamples, are strongly affected by the presence 
of the host; on the other hand, massive groups and groups with i? < do not show 
dependence on the host-group distance. Using a sample of groups identified in the final 
version of the 2dF Galaxy Redshift Survey, we find a very similar velocity dispersion 
behaviour in the observational data compared to results in the simulation. We also 
study the dependence of the groups velocity dispersion on the host masses in both, 
observations and simulation; finding that the larger is the host mass the higher is the 
effect on its vicinity. 
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1 INTRODUCTION 

In hierarchical models for galaxy formation, structures grow 
due to the accretion of smaller systems which formed at ear- 
lier times. According to this scheme, these systems should 
not only be in the cluster itself, but also in its vicinity, har- 
bouring an important number of groups or small clusters 
of galaxies. In this sense, clusters of galaxies represent high 
density enhancements in the distribution of galaxies. 

On the basis of this scenario, it is natural to ask whether 
physical properties of structures depend on its proximity to 
high density regions, which, as previously mentioned, may 
be well represented by rich clusters of galaxies; or, if on 
the contrary, they evolve independently on the presence of 
denser regions. 

Several works have focussed on the evolution of sub- 
structures inside massive systems in simulations or rich 
galaxy clusters (Tormen et al. 1998, Ghigna et al. 2000, 
Taylor & Babul 2003, De Lucia et al. 2003); nevertheless, 
not much attention has been paid to structures in the outer 
regions of these systems. Hence, it is not well known if some 
of the behaviours of substructures inside the clusters, could 
also be applicable to groups in their periphery. 



Einasto et al. (2003) have investigated the proper- 
ties of loose groups in the Las Campanas Redshift Survey 
(LCLGs), in the vicinity of rich clusters of galaxies, including 
clusters from the Abell and APM catalogues. X-ray clusters 
and also a sample of the richest groups from the LCLG cat- 
alogue itself. By using an additional sample, consisting of 
all those LCLGs which do not neighbour a rich cluster, they 
have been able to compare the properties obtained for dense- 
environment LCLGs with those of typical LCLGs. They 
found that in most cases, the observed richness of groups 
near rich clusters is larger than the corresponding value for 
groups in the comparison sample. The same effect is shown 
when using Abell counts as a measure of groups richness. 
The harmonic radius and the velocity dispersion are also 
somewhat larger in the neighbourhood of rich clusters than 
in typical loose groups. They also found a strong mass seg- 
regation; indicating that loose groups in the vicinity of clus- 
ters have masses that are larger than mean masses of groups 
drawn from the comparison sample. The same trend is ob- 
served in the luminosities of LCLGs, showing that groups 
around rich clusters seem to be more luminous. At last, they 
found that these properties of loose groups do not show any 
dependence on clusters richness. 
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Lemson & Kauffmann (1999), also address the issue ol 
properties of dark matter haloes as a function of local den- 
sity; concluding that mass is the only halo property which 
correlates significantly with local environment. In fact, high 
mass haloes arc undoroprcscntod in low density regions, 
whereas in high density regions the situation is the other 
way round. 

In this work we analyse the properties of groups in the 
vicinity of rich clusters; taking into account their masses, in- 
ternal energies, velocity dispersions and virial radii. To this 
end, we use data extracted from a high resolution numeri- 
cal simulation and from the final version of the 2dF Galaxy 
Redshift Survey (2dFGRS) (Colloss ot al. 2003). The outhnc 
of this paper is as follows. In section 2 we describe the simu- 
lation that wc used and give a brief description of the group 
identification algorithm; in section 3 we give a description 
of the 2dFGRS and the identification of groups in this cat- 
alogue; section 4 contains the results obtained for both, the 
simulation (subsection 4.1) and the observational data (sub- 
section 4.2); a summary and a discussion of our results are 
presented in section 5. 



2 NUMERICAL SIMULATIONS 

In this study we use one of the publicly available Very 
Large Simulations (VLS) (Yoshida et al. 2001). Given the 
success of the model to reproduce observational data, we 
have chosen the ACDM simulation, which has a box size 
of 479 Mpc h~^, and 512^ particles, resulting in a mass per 
particle 6.86 x 10^" Mq . The cosmological parameters 
of this simulation are: Qo = 0.3, Qa ~ 0.7, normalisation 
as — 0.9 and h = 0.7. Due to the large box size it is possi- 
ble to find a large number of rich systems, which give us the 
possibility to work with a reliable statistical sample. In addi- 
tion, the simulation resolution enable small systems to have 
enough number of particles to confidently estimate physical 
properties. 

The groups identification is carried out using a stan- 
dard friends-of-friends algorithm. The analysis of the VLS 

haloes catalogue is performed using a linking length given 
by / = 0.17n~^'''^, where n is the mean numerical density of 
particles, corresponding to an overdensity of 5p/p = 200. 
This criterion leads to ~ 379300 identified systems with 
masses ranging from 1.4 x 10^^ Mq (20 particles) to 
2.35 X 10^^ Mq h-^. 



3 OBSERVATIONAL DATA 

Observational clusters and groups used in this work arc ex- 
tracted from the final version of the 2dFGRS, with the best 
redshifts estimates within the Northern (—7°. 5 < 5 < 2°. 5, 
9''50'" < a < 14''50™), and southern (-37°.5 <5< -22°.5, 
21'' 40™ < a < 3''30'") strips of the catalogue. This galaxy 
sample contains 221414 galaxies. The group identification 
was performed using an algorithm similar to that developed 
by Huchra & Geller (1982), modified in order to take into ac- 
count sky coverage problems (Morchan & Zandivarez, 2002), 
such as redshift completeness, which represents a ratio of 
the number of galaxies for which redshifts have been ob- 
tained to the total number of objects contained in the par- 



ent catalogue; and magnitude limit mask, which corresponds 
to variations of the parent survey magnitude limit with the 
position on the sky. Also, it was taken into account a mi- 
nor efi^ect introduced by the magnitude completeness mask 
(/i-mask). Mask software was obtained from the final ver- 
sion data release of the 2dFGRS, kindly provided by Peder 
Norberg and Shaun Cole. 

In the construction of the group sample, values of 
5p/p = 80 and Vb = 200 km were used to maximise 
the groups accuracy, as stated by Mcrchan & Zandivarez 
(2002). The group catalogue comprises a total number of 
6088 galaxy groups with at least 4 members and mean radial 
velocities in the ranee 900 km s-^ < V < 75000 km s'^ . 
The minimal rmmber of members is imposed in order to 
avoid pseudo-groups. 

Virial masses were estimated using the virial radius and 
the velocity dispersion {Myir = cr^Rvir/G, Limber & Math- 
ews, 1960). For a better efficiency and stability in the estima- 
tion of the velocity dispersion, we use the biweight estimator 
for groups with more than 15 members, and the gapper esti- 
mator for smaller groups (Beers, Flynn and Gebhardt, 1990; 
Girardi et al., 1993; Girardi & Giuricin, 2000). Finally, the 
catalogue mean velocity dispersion is 253 km s^^ , the mean 
virial mass of 7.2 X 10^=^ Mq h-^ and a mean virial radius 
of 1.04 Mpc h-^. 



4 GROUPS IN THE VICINITY OF CLUSTERS 

In order to analyse the behaviour of systems surrounding 
high density regions, we discriminate "host" clusters from 
"groups" using a mass threshold Mcut- The studies were 
carried out in both, simulations and observational data de- 
scribed in previous sections. 

4.1 Groups properties in the simulation 

Host clusters in the VLS, were selected as systems with 
masses M > Mcut = 4 x 10^* Mq h~^, whereas systems 
with masses 1.4 x 10^^ Mq < M < Mcut are considered 
groups. With this criterion wc find ^ 379000 groups and 330 
hosts, leading to a mean host separation of ~ 70 Mpc . 
In this way we have defined our sample of groups neigh- 
bouring rich clusters, enabling us to study the dependence 
of groups properties on the host-group distance. 

The leading aim of our studies, is the dependence of 
main groups physical properties on the host distance. Given 
the asymmetry of these property distributions, we use the 
median instead of the standard mean. We find an increasing 
trend for the median velocity dispersion when considering 
smaller host-groups distances. Top panel of figure 1 shows 
the above mentioned trend, for all systems in our group sam- 
ple. As we are considering host-group distances starting at 
3 Mpc . miscalculation of group velocity dispersion ow- 
ing to hosts particles contamination, is not likely to happen. 
Assuring that the observed effect is not produced by high 
velocity host particles misplaced in nearby groups. 

As stated by Lemson & Kaufi^mann (1999), there is a 
dependence of the halo mass function on the environment, 
which is skewed towards high mass objects in overdense re- 
gions. Hence, the rising of the velocity dispersion in denser 
environments, could be the result of a higher abundance of 
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Figure 1. Host-groups distances against median velocity disper- 
sion. In the top panel the sample corresponds to all groups, in 
the middle panel to groups with masses 1.4 X 10^^ Mq h^^ < 
Mgroup < 2.0 X lO-"""' Mq h^^, and in the lower panel groups 
masses are 2.0 X lO^-^ Mq < Mgroup < 4.0 X lO" Mq h'^. 
Error bars were estimated by the bootstrap resampling technique. 



massive haloes. In order to test this possibility, we com- 
pute the mass function in 1 Mpc h^^ thick spherical shells, 
centered on hosts for several radii (3.5, 9.5, 20.5 and 40.5 
Mpc h~^). This allows us to quantify the dependence of the 
mass function on the overdensity (5= 5.89, 1.09, 0.28, 0.05 
respectively, for the above radii). Figure 2 shows in dashed 
lines the mass function for the mentioned radii multiplied 
by (1 + <5)~^ and labelled with the associated S's; solid line, 
represents the total mass function. Due to the overlapping, 
we shift the curves corresponding to the three highest over- 
densities by factors 10, 100 and 1000 respectively. As can be 
seen, the mass function shape does not change for groups 
with masses 1.4 x IO^^Mq hr'^ < M < 2 x lO^M© 
(vertical dashed lines in figure 2). Hence, the results ob- 
tained from this low mass sample will not be affected by the 
overabundance of massive haloes in high density regions. 

Based on the previous analysis, we explore the ob- 
served velocity dispersion behaviour for groups according 
to their masses; considering masses greater and lower than 
2.0 X 10^^ Mq h~^. As a result of this resampling we find 
that velocity dispersions of low mass groups are strongly af- 
fected by the presence of host clusters (middle panel of figure 
1), whereas high mass groups do not show a significant vari- 
ation of the median when approaching to their associated 
hosts (lower panel of figure 1). According to this, we focus 
the following analysis on this low mass group sample. 

At first sight, the rising of the median velocity disper- 
sion, could be a result of the gravitational influence of the 
host on its neighbourhood. This possibility encourages us 
to analyse the behaviour of the internal energy, E, in the 
same way we did for the velocity dispersion. This quantity 



Figure 2. Mass function around hosts, corresponding to dif- 
ferent overdensities 5. The solid line corresponds to the total 
mass function and dashed lines to the mass function computed 
in 1 Mpc h^^ thick spherical shells at radii 3.5, 9.5, 20.5 and 
40.5 Mpc h~^, multiplied by (1 + 5)~^ . For clarity, the three 
first curves and the total one, have been shifted (multiplied by 
1000,100 and 10 respectively) in the y-axe. Vertical dashed lines 
indicate the mass range of the low mass group sample. Error bars 
correspond to poisson estimation. The total mass function have 
error bars smaller than the point size. 



is estimated by: 



E = la^Mg 
2 



Gm'^N{N - 1) 



where the first and second terms correspond to the kinetic 
and potential energies respectively, a is the 3D velocity dis- 
persion, Mgroup the group mass, m the particle mass, N the 
number of members of a group and R^ir the group virial 
radius. Top panel of figure 3 shows the dependence of the 
group median energy on the host distance. It is clear that 
groups situated closer to high density regions, tend to have 
higher energies, indicating that these groups are typically 
weakly bounded. Next, we subdivide the low mass group 
subsample into positive and negative internal energy groups. 
In this case, we find that the rising median velocity disper- 
sion trend, is maintained for positive energies groups (open 
circles in bottom panel of figure 3), whereas the negative 
energies ones, show a constant behaviour (filled triangles in 
bottom panel of figure 3). Moreover, positive energy groups 
not only arc the main contributors to the signal, but also 
dominate in number towards central regions; as can be seen 
in the middle panel of the same figure. This result suggests 
that the rising of the median velocity dispersion towards 
smaller host-group distances, has to do with an equilibrium 
loss, due to an increment in the internal velocities of these 
groups. 

To complement the previously obtained results, we have 
also analysed the dependence of the median virial radius on 
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Figure 3. Lower panel shows host-groups distances against the 
median velocity dispersion. Dots: low masses group sample. Open 
circles: groups with energy E > 0. Triangles: groups with energy 
E < 0. Error bars were estimated by the bootstrap resampling 
technique. In the middle panel is plotted the ratio between the 
number of groups with E > and E < 0, as a function of the 
host-group distance. In this case error bars correspond to error 
propagation technique. In the top panel is shown the median en- 
ergy against host-group distance. Bootstrap resampling technique 
was used to obtain the error bars. 




Figure 4. Median virial radii (Top panel) and median mean den- 
sity (Bottom panel) against host-group distance. Error bars are 
estimated using the bootstrap resampling technique. 



Figure 5. Upper panel: median velocity dispersions for Mcut = 
4 X 10^^ Mq h-1,6 X 10" Mq /i-l and 9 x 10^* Mq h'^ cor- 
responding to the bottom, middle and upper trend respectively. 
Bottom panel: median velocity dispersion in the nearest bin of 
the host-group distances {(irnax), against Mcut- Numbers beside 
points indicate the number of host with masses larger than Mcut ■ 
Error bars are estimated according to the bootstrap resampling 
technique. 



the host distance. As shown in the top panel of figure 4, the 
median virial radius weaJtly increases towards inner regions. 
This is in agreement with the previous analysis, since con- 
tributes to a decrease of the potential energy with the conse- 
quent increment in the internal energy. In the bottom panel 
of figure 4, the median density of groups, p = (N/Rl^^), is 
plotted against the host-group distance showing that groups 
nearer to hosts, have lower mean densities than groups in the 
outer regions. 

A possible interpretation of the above results can be 
found analysing the dynamical properties of systems of 
galaxies in the presence of external gravitational fields. The 
interaction between two systems of very different size will 
mainly affect the properties of the smaller one. The tidal ra- 
dius Rt specifies the size of the smaller system beyond which 
the mass becomes bound to the potential of the larger one. 
For a group with total mass Mg in the neighbourhood of a 
host with mass Mh, the tidal radius can be roughly approx- 
imated by 



Rt 



\3MhJ 



where r is the host-group distance (Binney & Tremaine, 
1987). For our low mass group sample, the mean group mass 
is ~ 4.0 X 10^^ Mq h~^; whereas host sample has a mean 
mass of 6.5 x 10^'' Mq h^^ . Replacing these values in 
the above equation, for a distance of 3 Mpc , a, value of 
Rt = 0.38 Mpc is obtained, which is quite similar to 
the median virial radius at the same distance (top panel of 
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figure 4). This result siiows tiiat in the proximity of large 
systems, tidal forces exerted by the host group can be as 
strong as groups internal forces. 

In an attempt to explore the dependence of groups 
properties on the hosts masses, we compute median veloc- 
ity dispersions in the same way we did before, for differ- 
ent mass thresholds. In this case we leave a gap between 
the upper mass limit used to select groups, and the mass 
threshold used to select hosts, in order to work with the 
same low mass group sample we did before. In top panel 
of figure 5, we show the median velocity dispersions against 
host-group distance, for three different host samples using 
Mcut = 4xl0",6xl0^'' and 9x10^"* Mq corresponding 
to the lower, middle and upper trend respectively; the more 
massive are the selected hosts, the steeper is the median ve- 
locity dispersion trend. In the bottom panel of figure 5 we 
plot the value of the maximum median velocity dispersion, 
(Jmax, corresponding to groups in the distance range from 3.0 
to 3.45 Mpc , as a function of Mcut- A strong enhance- 
ment of the effect that hosts exert on its neighbourhood, 
is observed when considering more massive hosts. This fact 
supports the scenario where the gravitational potential of 
massive systems, perturbs the behaviour of the surrounding 
smaller groups. 

4.2 Groups properties in the 2dFGGC 

In this section we describe the results obtained for the sam- 
ple of groups derived from the 2dFGRS as described in sec- 
tion 3. 

In order to compare with the results of the simulation, 
we extract a random subsample from groups in the VLS, 
that reproduces the 2dF group velocity dispersion distribu- 
tion. Inset plot of figure 6, shows in dashed line the nor- 
malised distribution imposed for VLS groups, and in solid 
line the corresponding to 2dF groups. We use for both sam- 
ples Mcut = 7 X 10^'' Mq h~^\ this choice leads to 45 hosts 
and 4548 groups in the 2dF, and 156 hosts and 38034 groups 
in the VLS. We compute, as in previous section, the median 
velocity dispersion of groups; in figure 6 solid line repre- 
sents the median velocity dispersion for 2dF groups, whereas 
dashed line is associated to VLS groups. This figure shows 
an enhancement in the velocity dispersion, as predicted in 
the previous section using the total VLS group sample. 
The mentioned behaviour, is consistent with Einasto et al. 
(2003), since they state, that the velocity dispersion of loose 
groups in the neighbourhood of rich clusters is about 1.3 
times larger than that in field groups; value which is very 
similar to the ratio between the first and last point of figure 
6. They also find that groups virial masses in high density 
environments are about 2.5 times larger than masses of field 
groups, which is also in agreement with our results, since 
the virial radii for 2dF groups show a very weak increment 
towards smaller host-group distances, and the virial mass 
is Muir oc RyirfJ^. This results in an evident increment on 
the groups virial masses towards the hosts, mainly due to 
the velocity dispersion behaviour. On the other hand, we 
have seen in section 4.1, that real masses, do not have a de- 
pendence on the host-groups distance; indicating that virial 
masses might not be a reliable indicator of the real mass of 
a group nearby to a high density region. In order to study 
the dependence of groups velocity dispersion on the host 
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Figure 6. In the inset panel, the solid line represents the 
2dFGGC normalised velocity dispersion distribution; the dashed 
line is for a subsample of VLS groups which imitates the former. 
Velocity dispersion medians for both samples, are plotted in the 
main panel for Mcut = 7 X lO'^* Mq h~^, the line types cor- 
respond to the same samples as for the inset panel. Bootstrap 
resampling technique is used for the computation of error bars. 



masses for the 2dF sample, we take a higher mass thresh- 
old, Mcut = 9 X lO^"* Mq h~^, resulting in 25 hosts, and 
perform the same kind of analysis. Figure 7 shows a higher 
signal for the velocity dispersion in the inner regions which 
is in concordance with the results we found for the simula- 
tion; the larger is the host mass, the higher is the typical 
velocity dispersion of groups near to the host. 



5 CONCLUSIONS 

Given the large volume and high resolution of the publicly 
available VLS, we could identify a large number of systems 
of galaxies in a wide range of masses. Using a mass threshold 
we defined several samples of host systems and groups. 

Our results show that the median velocity dispersion 
of groups of galaxies increases as the distance to the host 
clusters is smaller. Taking subsamples of groups, consider- 
ing high and low masses, we find that the previous behaviour 
is not present in the high mass subsample whereas low mass 
groups account for the observed effect. This low mass sub- 
sample was used to study the behaviour of groups according 
to their internal energies. We found that the internal en- 
ergy of groups increases towards smaller values of host-group 
distances, showing up that the gravitational equilibrium of 
groups is infiuenced by the presence of massive systems. By 
splitting this sample in two groups subsamples with posi- 
tive and negative energies, we found that only those with 
positive energies increase the median velocity dispersion to- 
wards smaller host distances. In order to unmask the source 
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Figure 7. Median velocity dispersion for Mcut = 9 X 
10^* Mq Line types correspond to the same samples as for 
figure 6. Bootstrap resampling technique is used for the compu- 
tation of error bars. 



of these effects, we studied the dependence on the distance of 
complementary physical properties. By considering smaller 
host-groups distances, the virial radius increases and the op- 
posite trend is observed for the mean density. 

A possible scenario describing the above results, may 
be placed within the framework of the tidal forces exerted 
by the massive clusters over neighbouring systems. In the 
proximity of a cluster, its potential has a strength which is 
similar to the one exerted by the group itself; indicating that 
tidal forces are likely to affect groups properties. 

Finally, analysing different subsamples of host clusters 
we found that, the more massive is the host cluster the 
greater is the effect on the group velocity dispersion, which 
is also consistent with the previous model. 

Using the final version of the 2dF redshift survey, we 
repeated the velocity dispersion analysis performed to the 
VLS groups. In this case we take a subsample of VLS 
groups, which reproduces the velocity dispersion distribu- 
tion of groups identified in the 2dFGRS. The trend we find 
for the observational data is in excellent agreement with the 
result we obtain for the simulation subsample. Moreover, 
2dF groups velocity dispersion show a dependence on the 
host masses very similar to that found in the simulation. 
The good agreement between the results coming from the 
numerical simulation and the observational data suggests 
that no astrophysical mechanisms, other than the gravita- 
tional forces, are needed in order to explain the dynamical 
properties of groups of galaxies. 
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